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Abstract. We consider equilibrium level occupation numbers in a Fermi gas 
with a fixed number of particles, ra, and finite level spacing. Using the method of 
generating functions and the cumulant expansion we derive a recurrence relation 
for canonical partition function and an explicit formula for occupation numbers 
in terms of single-particle partition function at n different temperatures. We 
apply this result to a model with equidistant non-degenerate spectrum and obtain 
close-form expressions in terms of g-polynomials and Rogers-Ramanujan partial 
theta function. Deviations from the standard Fermi-Dirac distribution can be 
interpreted in terms of a gap in the chemical potential between the particle and 
the hole excitations with additional correlations at temperatures comparable to 
the level spacing. 
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1. Introduction 

Applications of statistical mechanics to fermion systems with discrete spectrum, such 
as semiconductor quantum dots [U [2], naturally involve single-particle averages in 
statistical ensembles with a fixed number of particles. In particular, the kinetic 
theory of tunneling [3j |4] through quantum dots with fast intra-dot electron relaxation 
involves average level occupation number, 

pZn dek 

in a Gibbs distribution of n independent fermions populating a set of single-particle 
energy levels {tk} (enumerated by fc = 0, 1, 2 . . .). Here Z„ is the canonical partition 
function 

Zn = ^ exp^-/3 ^ Vkekj ^n, J]^ , (2) 

and P is the inverse thermodynamic temperature. For fermions, occupation numbers 

Vk in the sum ^ take values and 1. 

The behavior of {vk)n is simple in two extreme limits of the typical level A. For 
<C 1 and large n, Eqs. ([I])-© reduce to the standard Fermi-Dirac distribution, 

{vk)n -> f{ek - A*) = ^ ^ ■ (3) 

Here (x is the chemical potential determined by the normalization condition n = 
/(sfc ~ m)- In the low temperature limit, /3A ^ 1, most of the statistical weight in 
Eq. ^ is in the ground state (defined as z/^ = 1 for < fc < n and Vk = otherwise). 
In this case it is common [B] to take only one excited state into account resulting 
in a two-state Gibbs distribution which is equivalent to Eq. for fc = n — l,n with 
/X = (e„ + e„_i)/2 and /3 -> /3* = 2/3 [3]. 

For finite /3A, the average occupation number {i'k)n deviates from Eq. in a 
non- universal way which depends on the details of the energy spectrum [21 [3 IB] ■ Exact 
analytical investigation of this regime is complicated by the combinatorial explosion in 
the number of levels with comparable statistical weights. Particle number projection 
technique based on Fourier extraction [91 [7] gives an exact closed-form formula for 
(j^fc)n (see Eq. (140) of [2 ) that scales quadratically with the number of levels; however, 
its potential for analytical investigation appears to be limited due to the additional 
sum over the Fourier variable. 

In the paper we address the problem of exact evaluation of canonical occupation 
numbers {i'k)n by deriving a general formula that scales linearly in the number of 
levels and quadratically in the number of particles, Eqs. (|9]) and (jTTI) below. The 
energy spectrum {e^} enters the formula only in terms of the single-particle partition 
function Zi computed at n different temperatures. We apply this general result to 
equidistant spectrum, ek = fcA, and derive an exact formula for {i'k)n in terms of 
polynomials in q = e~^'^, Eq. p^ . In the limit of degenerate Fermi gas, n[3A ^ 1, 
these polynomials converge to partial theta function [T0\ which is involved in a number 
of combinatorial proofs [TTJ[T^[T3] of Ramanujan's identities [T?]. The exact result for 
the equidistant spectrum can be approximated well by tailoring two standard Fermi- 
Dirac distributions ([3]) with different chemical potentials for holes and for particles, 
IJ'h = £n and fip = e„_i, respectively. In the high- and the low-temperature limits, 
this approximation converges to the asymptotically exact Fermi-Dirac and two-state 
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Gibbs distributions, respectively. At intermediate temperatures, /3A ~ 1, particle-hole 
correlation effects due to fixed n result in finite deviations from the exact solution. 



2. General expressions for fermion partition functions and the occupation 
numbers 

Grand canonical partition function Y serves as a generating function for the canonical 
partition functions Zn if expanded power series of z = e^^° , 

oo 

^exp(-/3^z.fe(efc-Aio)) =1 + J2z^z". (4) 

{i/fc} k n=l 

Y{z) is most conveniently calculated via its logarithm |15j . 

oo 

lny(2)= ^ln(l + ze-''^'=) -^^z", (5) 

k n—1 

where k„ ee (-l)"+i(n - l)!Zi(/3n), and 

k 

is the canonical partition function of a single particle. 

Relation between nlZn and k„ is the same as between the raw moments and the 
cumulants of a univariate probability distribution and is given by the complete Bell 
polynomials [16 , 

Zn = iniy^B{Ki,K2,...,Kn)- (7) 

The latter satisfy a recurrence relation [17], 

/n - 1\ 

), (8) 



m— 1 

which translates into 



n 

Zn = -y] (-l)™+iZi(/3m)Z„_™ . (9) 



n 

m— 1 



We set Zq — 1 identically. 

Combining Eqs. ([T]), (|4]) and ([5]) gives the generating function for the occupation 
numbers: 

^W,Uz^z- = f . (10) 

n—l 

Expanding the r.h.s. in powers series in z gives 



1 " 

Wk)n = ^ E . (11) 

m— 1 

Equations ([5]) and (fTTj) constitute our main general result. 
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3. Example: equidistant spectrum 

3.1. Exact finite-n results: polynomials 

For efc = /c A, the grand canonical partition function ([S]) can be expressed by an infinite 
product, 

oo 

Y{z) = l[il + q''z) = {-z:q)^, (12) 

k=0 

where q = e~^^ and (■;'?)« is the (/-shifted factorial [l8l 119) . 

Using g-analog binomial theorem of Euler (T^, formula 17.2.35) we can get the 
partition function directly from Eq. 

Zn^ \ , ■ (13) 

Applying Eq. ([TT]) . and transforming (/-shifted factorials ([T^, formula 17.2.13), 
one gets 

{iyk)n = l-p{k,n;q) , (14) 

n 

p(fc,n;(/)^ ^(7"('=+i) ((/-"; (Z)™ (15) 

n m — 1 

= 1 + ^ [| (^fe+l _ ql+k^n+l-^ (^g) 
m=l /=0 

Equation ()14|) defines occupation numbers for n fermions populating equidistant levels 
at equilibrium. It is clear from the explicit form (|16p that p{k,n;q) is a Laurent 
polynomial (the product contains negative powers of q if k < n). However, since 
< (i'fc)n < 1 for g by definition ([1]), the negative powers of q must cancel, thus 
we conclude that p{k,n;q) is always an ordinary polynomial in q for n > 0, fc > 0. 
This cancelation is not trivial and deems further mathematical investigation \20\ . 

A number of recurrence formulas can be derived for p(k,n\q) [3T], including a 
symmetry relation 

q''p{k,n;q) = q^p{n,k;q) . (17) 

Using (IT51) in the r.h.s. of ([T71) gives a sum of products with no negative powers of q 
at k < n. 

3.2. Largc-n limit: q-analog of Fermi- Dirac distribution 

If e„ P""^ then the Fermi gas is degenerate [TS] and the limit of (7" is 
appropriate. For k,n ^ 00, and k — n = const > —1 the polynomial sum in Eq. (|16p 
becomes a geometric series which gives 

jim p{k, n; q) - 0(-(?^-"+l/^ q^/^) , (18) 

k>n-l 

where 

CxD 

0((j,,7)= ^a™,/"' (19) 

m=0 
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Figure 1. Comparison of different approximations for the canonical occupation 
numbers in a Fermi gas witfi equidistant spectrum, n = 4 and /3A = 1. Crosses 
mark tfie exact values, tfiick (blue) continuous line shows the q-analog of Fermi- 
Dirac distribution, Eq. I|20|l . thin continuous line — substitution /q — > / in 
Eq. ||20|I . and short-dashed line — single Fermi-Dirac distribution, Eq. ((3)l. The 
long-dashed lines between and e„ are extrapolations of the corresponding 

functions into the gap /ip < e < /i/i. 



is known as partial theta functioifl [TU]. The partial theta function is famous for 
a number of identities discovered by Ramanujan in his lost notebook [M]. These 
identities have been studied extensively [22l [10], including some recent proofs by 
combinatorial methods [HI [121 [l3] . Using the symmetry relation (fTT]) for k < n, gives 
the particle-hole complementary result : Muin-^oo pik,n; q) = 1 — q"^'^^^/^, g^/^). 

In terms of level energies e^, the occupation numbers in a canonical degenerate 
Fermi gas with constant levels spacing A can be written as 

where 

f,{e)^9{-e-P\q'/^)'^' f{e). (21) 

and fj, = {en + e„_i)/2.The function defined in Eq. ([20| can be considered a g-analog 
[18] of the standard Fermi-Dirac distribution ([3]) since lim^^i fq{e) = /(e). 

Equation (|20l) expresses two essential deviations of canonical occupation numbers 
from Fermi-Dirac distribution. Firstly, the q- analog fq is different from /. An 
approximation of substituting fq—^f in Eq. ([20]) becomes exact both for g — > 1 and 
for g— >0. Numerically, we find maximal absolute deviation |/(efe) — fq{ek)\ of 0.0567 
reached for k = n, n—1 at /3A = 0.752. For the two levels closest to the gap, k ~ n—l, n 
approximating fq^f gives {vn)n = 1 — {vn-i)n = exp(— /3A/2)/[2 cosh(/3A/2)] which 
is equivalent to two-state Gibbs approximation [3]. A comparison between the exact 

t Note that lim„-^oo n; q) = Y.m=0^~'^)"' l"^'""^'^^^ ^'^ instance oi false theta series [22] in 

the sense of L .J. Rogrers [23) . 
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result ([H]), the large-n limit ([20]) and a single Fermi-Dirac distribution is shown in 
Fig. [U for n = 4 and ;3A = 1. 

Secondly, barring the difference between fq and /, Eq. (PH)) can be seen as a 
combination of two Fermi-Dirac distributions with different chemical potentials for 
particles, /Xp = /x — A/2 ~ e„_i, and for holes, — /i + A/2 = e„, respectively. This 
is precisely what is to be expected if one considers particle and hole excitations from 
the ground state as two uncorrelated microcanonical ensembles. In this view, the 
moderate difference in the functional dependence between /(e) and fq{e) results from 
correlation between a particle and a hole created in a single pair-excitation act. Note 
that the difference //ft, — /Xp = A can not be ignored even as /? ^> oo and g — > 0, thus 
using a single Fermi-Dirac distribution, Eq. ([3]), necessarily fails for finite A and low 
temperature. 
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